Holding strategies in a bus-route model 



Scott A. Hill 
Department of Physics, Southern Methodist University 
(Dated: February 5, 2008) 

A major source of delays in public transportation is the clustering instability, which causes late 
buses to become progressively later while the buses trailing it become progressively earlier. In this 
paper, we study this instability and how to neutralize it using the common practices of holding 
and schedule slack. Starting with an on-time route, we delay one or more buses at a single stop, 
and determine how these delays grow over time. We compare the effects of two different types of 
holding on the stability of the system, and briefly investigate how our results change with the use 
of timepoints. 



I. INTRODUCTION 

The typical bus or train system is intended to provide 
regular, periodic service, so that each stop has an ar- 
rival every n minutes. In practice, however, public trans- 
portation often suffers from maddening inconsistencies, 
which discourages potential passengers from depending 
on it for their daily commute. Given the environmen- 
tal, economic, and political impact of increased gasoline 
consumption, it is important for us to understand these 
delays and reduce their frequency, if possible. 

One cause of these delays is an inherent inst ability pQ 
in the dynamics of a bus route: when a single bus (bus 
A) is delayed, it has to pick up more passengers which 
delays it further, while the bus following it (bus B) has 
fewer passengers to pick up, and runs faster. In many 
cases, bus A is unable to recover and, on busy routes, 
bus B may even catch up to bus A, forming a cluster. 
This is not always a problem: for an evening outbound 
route, where most passengers board at the beginning of 
the route, it does not matter that the service is irregular 
so long as the passengers are gotten to their destinations 
quickly. However, in most circumstances, this instability 
results in longer, unpredictable passenger waits. 

Bus dispatchers can counter this behavior with holding 
and slack. Buses which are running too fast (such as bus 
B above) are held at a stop so that they do not catch up 
to the bus in front. In schedule-based holding, buses are 
prevented from leaving a stop until a scheduled departure 
time, while in headway-based holding, buses are prevented 
from leaving until the preceding bus is far enough away. 
Most bus services implement holding only at a few time 
points along the route, while light rail services, which 
typically stop at every station anyway, can implement 
holding at every stop. 

If holding is only implemented when buses are running 
more quickly than usual, then it is only partially effective 
in preventing clusters: it may keep bus B from catching 
up to bus A, but it does not allow bus A to recover from 
its delay. However, if most buses are held, then not hold- 
ing a delayed bus may be enough to allow it to recover. 
This is done by introducing slack into the schedule; that 
is, by allowing more time for buses to travel from stop to 
stop than is ordinarily necessary, so that the typical bus 



is held at every stop. 

The clustering instability was first mentioned in the 
engineering literature by Welding [1 in 1957; since 
then, it has been covered extensively in the engineering 
literature [21 U 0], although not to the point of exhaus- 
tion. A few physicists, inspired perhaps by the much 
larger field of traffic study [5], have turned their attention 
to modelling bus routes [6j [7j |9j [10], using a number 
of methods (cellular automata, linear stability analysis, 
and chaos theory). 

The engineering literature has primarily focused on the 
statistical properties of multiple delays along an entire 
bus route. In this paper we take a more "microscopic" 
approach, by investigating the deterministic propagation 
of delays which occur at a single stop, on a bus route 
that is otherwise running on time. Using a model which 
is similar to others in the literature [U HI EH] , we calculate 
the maximum delay that a bus can recover from (a quan- 
tity we call the "buffer") for a given amount of slack, 
and how this buffer varies when two or more consecu- 
tive buses are delayed at the same stop. We show that 
headway-based holding is typically better in recovering 
from random delays, although schedule-based holding al- 
lows buses to recover more quickly when the delays are 
shorter. We mostly assume that buses are held at every 
stop; however, we end by briefly showing the effect that 
timepoints have on the buffer of a single bus. 



II. MODEL 
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FIG. 1: How buses and stops are indexed in this model. 



Consider a series of irregularly spaced bus stops, la- 



2 



belled with the index s, and a series of buses labelled 
with the index b. Each bus visits each stop in increasing 
order, and each stop is visited by each bus in increasing 
order (Fig. [T]) . We define tb, s to be the time (in minutes) 
at which bus b departs stop s. Before bus b can depart 
stop s, it must do three things: 

1. depart stop s — 1, at time i; 

2. drive from stop s — 1 to stop s, which takes time T s 
; and 

3. pick up passengers at stop s. 

As a simplification, we ignore the time it takes for passen- 
gers to get off the bus: this corresponds to the situation 
in which the number of alighting passengers is either very 
small, or in which passengers can alight through a rear 
door at the same time other passengers are boarding. 

The time it takes to pick up passengers at stop s is 
equal to the time it takes each passenger to board (the 
unit boarding time), multiplied by the number of passen- 
gers waiting at the stop. This number of passengers is, 
in turn, equal to the wait since the last bus (tb jJS — i,s) 
divided by the time it takes each passenger to arrive (the 
interarrival time, the reciprocal of the arrival frequency) . 
Thus, the time it takes to board passengers is 

boarding time = /i(tf,. s — tf,_i iS ), (1) 

with the passenger constant \i > defined as 

unit boarding time 



I 1 = 



(2) 



interarrival time 
We can now specify a recursion relation for i& iS : 

tb, s = h,s-i + T S + fi(tb, s — i&-i,s)- (3) 
This can be simplified with a change of variable 

s 

h.s —> h,s + ^""^ Tj, which eliminates T s from the equa- 

i=0 

tion entirely. If we then solve Eq. [3] for tb, s (and T s re- 
moved), we have 



tb,a = (1 + v')tb, s -i 
where we define 



n't b - 



1-M 



(4) 



(5) 



To implement holding, we must define a schedule func- 
tion Sb, s , which specifies the time at which bus b should 
depart from stop s. For simplicity, we assume a homoge- 
neous schedule in which buses are evenly spaced: 



S b , s = bA + s(fiA + a). 



(0) 



The parameter A is the time between successive buses 
at a particular stop, fj,A is the time it takes a particular 
bus to travel from stop to stop under normal conditions, 



and a is the slack built into the schedule [TT] • If there is 
no slack (<r = 0), then tb, s = Sb, s is a solution to Eq. El 

We implement holding by specifying some earliest de- 
parture time i™ n ; if the bus is ready to leave before i min , 
it is held so that it leaves at i mm . We do this by rewriting 
Eq. [4] as the conditional 



(1 + n')t b „ 

^min 

l b,s 



n'h- 



(7) 



Note that our model will hold buses at every stop, rather 
than at the less frequent time points seen in most bus 
systems. 

It is convenient to work with the delay of a bus at any 
given stop, rather than its departure time. We define the 
unnormalized delay of bus b at stop s to be 

h,s = h,s — Sb, s , (8) 
and the normalized delay (or, simply, the delay) to be 
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(9) 



Delay is measured relative to the schedule function Eq. [6] 
early buses have d < 0. We rewrite Eq. [7] in terms of 
delay, to arrive at our main dynamic equation 

d b , 3 - max | { ^' )dh ' S ~ 1 ~ fX ' db ~ 1 ' S ~ ^ ( 10 ) 

where cdb-\. s is the minimum delay of bus b at stop s. If 
c = 0, then buses are never allowed to be early (and the 
delay db, s is never negative), and we have schedule-based 
holding. If c = 1, then buses are never allowed to be 
earlier than the bus preceding them [db s > c?b_i jS ), and 
we have headway-based holding. When the second of the 



two conditions in Eq. 10 is larger, we say that "holding 
has been triggered" . 

In the calculations that follow, it will be necessary to 
solve several recursion relations of the form 



x s — zs s _i + Az s + 7; 
this has the solution 
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III. SOLUTIONS 

We now consider a case in which, following a series of 
on-time buses (db ;S — for s < and b < 1), one or more 
buses, starting with bus b — 1, are delayed at stop s = 0. 



The First Bus 



Given our assumption that db 



for b < 1, the 



dynamic equation (Eq. 10 1 for the first bus is 
d M+ i = max [(1 + fi')di, s - fi', 0] . 



(13) 



3 



If the bus is early or on-time at any stop so (di iSo < 0), 
it will be on-time from then on (c?i jS>So =0). As long as 
d± tS > 0, however, Eq. 13 has the solution (cf Eq. 12 1 

di, s = 1 - (1 + M') S (! - di,o) (14) 

which grows exponentially towards positive or negative 
infinity, depending on the sign of 1 — d\a. Thus there are 
two types of results (as shown in Fig.|2|: 

(a) If d\xi < 1, then the delay of the bus decreases until 
it is on time (di. So < 0), after which the bus remains 
on time. We say that such a bus is recovering. The 
bus has fully recovered when 



_ -ln(l-rf 1|0 ) 
S>S0= ln(l + M ') • 



(15) 



Note that bus routes with larger numbers of pas- 
sengers (larger values of //) recover more quickly 
(Fig. |2|. A special case of this is when sq < 1, 
in which case the bus recovers instantly; this occurs 
when dtp < //. 

(b) If difl > 1, then the delay of the bus increases ex- 
ponentially; we say that such a bus is unrecoverable. 
A single unrecoverable bus delays all following buses: 
either the bus behind it catches up to it (if c = 0), 
forming a cluster, or else it also becomes exponen- 
tially late (if c — 1); in either case, the homogeneous 
behavior has broken down irrevocably. 

We define the buffer (3 of a bus to be the largest delay 
df,fl that it is able to recover from; the buffer of the first 
bus is thus Pi — 1. This is a normalized quantity based 
on the normalized delay db, s - In terms of the unnormal- 
ized delay £b,s (Eq. pi, the bus recovers instantly when 
^b.s < cr, because the bus need only eliminate part or 
all of its slack. The bus's unnormalized buffer is a/ fi: as 
Newell [3 J explains it, if the slack per stop a is larger than 
the time [ilh,s it takes to board the additional passengers 
due to the delay, then the bus will eventually recover. 



B. The Second Bus 

If the first bus is late, it will pick up some of the pas- 
sengers originally destined for the second bus, and so the 
second bus will run faster. Therefore, if the second bus 
starts out on-time, it will have no trouble remaining on- 
time, so long as the first bus is not unrecoverable. It also 
stands to reason that if the second bus is delayed at stop 
s = 0, then it will recover more quickly, and have a larger 
buffer (3, than if the first bus were not also running late. 

The dynamic equation (Eq. 10 1 for the second bus is 



(1 + /j,')d 2<s -i - fjf(d 1>a 
cd\, s 



1) 



(16) 



If the first bus is unrecoverable, then the second bus will 
be unrecoverable as well, so we need only consider the 
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FI G. 2 : The normalized delay of the first bus, calculated from 
Eq. |l3| In Fig.Sk, the passenger constant is fixed at // = 0.1, 
while the initial delay is varied: examples of unrecoverable, re- 
covering, and instantly recovering buses are shown. In Fig. ^p, 
the initial delay is fixed at di.o — 0.9, while the passenger con- 
stant is varied. Notice that the bus recovers faster when the 
passenger constant is greater. 



case where di.o < 1. If d\$ < /i, then the first bus recov- 
ers instantly, and the behavior of the second bus can be 
described using Section flll A| Thus we only need consider 
the case where < di j0 < 1 (note that = fx). 

There are two cases to be considered: 



(a) If the holding condition in Eq. 16 is triggered at some 
stop s— 1, then d2,s-i — cdi,s-i- Now as long as the 
first bus hasn't recovered, Eq. [14] can be rewritten as 



rfi. 



and so 

d 2 .s 



(d liS +//)/(! + //) (17) 
/i'di, s — n(l — c) 



cdi 
cdi 



(18) 



Clearly the upper term is smaller than the lower term, 
and so the holding condition is triggered at stop s as 
well. Therefore, once the holding condition has been 
triggered for bus 2, it will continue to be triggered 
from there on out. If c = 0, then the second bus 
recovers immediately; if c — 1, then the second bus 



4 



2.2 



1.6 



1.4 



1.2 - 



X 1 



■~ 0.8 - 

N 

"I 0.6 - 

F 


z 0.4 



































Bus 
Bus 2 


1 recovers 
unrecoverable 




jf 














































!<* 

I.OOl 

0.1 ■■ 






'us 1 
overabl 












l 

unrec 


e 








0.2 -- 
0.5 - 












Both buses recover 











































0.2 0.4 0.6 0.8 1 1.2 1.4 
Normalized initial delay of first bus d, 

FIG. 3: A diagram showing the buffer of the second bus as a 
function of the first bus's initial delay di,o. One can interpret 
this as a phase diagram: in the unshaded region, both buses 
recover, while in the shaded regions one or both buses are un- 
recoverable. The buffer depends on the passenger constant fi' , 
being smaller for larger passenger constants, but never going 
below 1. This figure is independent of the holding strategy c. 



recovers along with the first bus. In either case, how- 
ever, once the holding condition is triggered, both 
buses recover. 



(b) If the second bus's holding condition is never trig- 
gered, and the first bus hasn't recovered yet (i.e 
s < s ), we can substitute Eq 
get 



14 



into Eq. 16 



to 



d 2 .s - (1 + fJ.')d 2 , s -i + + M') S (1 ~ di,o) ~ 2^', (19) 
which, according to Eq. [12] has the solution 
<h,. = 2 + [//(l - d lfi )s - (2 - da,o)] (1 + /"T- (20) 

Once s > so, di l8 = and the second bus follows 
the "first bus" pattern; therefore, the second bus will 
ultimately recover only if d2,s < !• Since 



l 2,S 



// ln(l - difl) 2-d 2 .o 



ln(l + //) l-di,o' 
the inequality e?2.s < 1 holds if <i2,o < (3% where 



(21) 



#2 = 1 + di i0 + 



ln(l + /x') 



-(l-d li0 )ln(l-d lj0 ). (22) 



Figure [3] shows a graph of the buffer of the second bus 
as a function of the first bus's initial delay. The buffer's 
curve marks the transition between the case in which 
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FIG. 4: Figure [4^, shows the buffer of the third bus as a 
function of the delay of the second bus, and for two values 
of the delay of the first bus (ti,o = 0.4 and ii,o = 0.8); all 
three buses recover in the region underneath the correspond- 
ing buffer curve. Clearly the buffer of the third bus depends 
on the holding strategy c. Figure [4]d shows how the three 
buses behave at the X marked in (a), with ti t o = 0.8. For 
schedule-based holding (c — 0), the third bus is unrecover- 
able, while for headway-based holding (c = 1), all three buses 
recover. The passenger constant is // = 0.1 throughout. 



both buses recover (below) , and the phase in which one or 
both buses are unrecoverable (above and to the right): we 
see that when the first bus is very late (but recovering), 
the second bus can recover from a much larger delay. 
As the passenger constant // increases, the region where 
both buses recover becomes smaller; this is because, for 
larger passenger constants, the first bus recovers more 
quickly (as seen in Fig. [2]) and so the second bus does 
not get as much benefit from the first bus's delay. 



C. The Third Bus 

Although the recovery of the first two buses is inde- 
pendent of the holding strategy c, the same is not true 
for subsequent delayed buses. Fig. [4ji-b shows the buffer 
of the third bus (6 = 3) as it depends on the initial delay 
of the first two buses; this clearly depends on the hold- 
ing strategy. Fig. Hfc shows why this difference exists: if 
schedule-based holding is in place, then the second bus 
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FIG. 5: The solution of Eg. [10] when db,o = T and c = 0; buses 
alternate above and below the steady-state solution, shown as 
a solid black line here. 



recovers quickly from a small initial delay, and the third 
bus loses the benefit of following a late bus. Note how, 
when c = 0, the third bus's buffer approaches 1 (the 
buffer of a bus following an undclayed bus) when t^o ap- 
proaches zero. When c = 1 and d 2j o is small and di<o is 
large, on the other hand, the second bus can only recover 
as quickly as the first bus, giving the third bus a larger 
buffer. In a sense, the delay of a recovering bus is a re- 
source which makes the trailing buses more resistant to 
delay. 



Then Eq. 10 becomes 

(1 + //)di )S _i - n'di, s - fj! 



dh.s — 



(23) 



If di_ s = and di )S _i = 0, then db, a = and 
the hypothesis is satisfied. If e?i jS _i but 
dx <s — 0, then (1 + (i')di,„-i < fJ,' (according 
to Eq. 13 1, which means that the first case in 
Eq. 23 is negative, and so db >s = = d± tS and 
the hypothesis is again satisfied. If d 1-s ^ 
and di, s _i ^ then, according to Eq. 
(1 + //)di jS _i = (ii. s + fi', and so 
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di, s +// - ^'dx >s 
di, a 



(24) 



and since (1 — fj/)di jS < c?i. s , the second con- 
dition applies and <4 s = d\ s . Q.E.D. 



The steady-state solutions are thus 



lim dh 







1+/*' " ~ " . (25) 

\-{1 + ^) s {1-t) c=l ^ ; 



It can be shown that, as long as d^i > for both cases, 
the c = case in Eq. [25] reaches zero faster than the 
c = 1 case: this means that schedule-based holding allows 
the buses to recover more quickly than headway-based 
holding. 



D. Many Buses 

To expand our model to the many-bus limit in a man- 
ageable way, we consider the case where all buses are de- 
layed by the same amount r at stop s — (i.e. dbo = t 
for all b > 1); this might be caused by construction or 
traffic on one particular street. If r > f3i = 1, the first 
bus, and eventually all buses, will be unrecoverable; thus 
we only consider the case when r < 1. Any bus with 
(4,0 < 1 will ultimately recover, so our concern is not 
whether this situation recovers, but how quickly it recov- 
ers. 

For c = (Fig. [5]), the solutions dt,. s approach (in an 
alternating manner) the steady-state normalized solution 
db, s = t — sfj, as b increases. This solution can be better 
understood by looking at the corresponding unnormal- 
ized solution ij, s — -r — sa: buses make up time by 
eliminating slack from the schedule, becoming a minutes 
earlier at each stop until they have recovered. 

For c = 1, we can prove that (4 iS = di. s for all s: 

Proof: by induction over b and s. The base 
case b = 1 is automatic, while the base case 
s = is true because d^^ = d i a = t. Sup- 
pose that db', s ' = di tS for all b' < b and s' < s. 



IV. TIMEPOINTS 

So far we have considered the situation where holding 
occurs at every stop. This may be a realistic description 
of a light-rail system, where trains stop at every stop 
anyway. However, a typical bus route may have a stop 
at every intersection, and a model which forces buses to 
stop and wait at each of these is unrealistic. In reality, 
most bus systems designate a few stops as timepoints, 
and only implement holding there. To model this, we 
assume that every TVth stop is a timepoint, and evaluate 
Eq. [10] for 6=1, with the stipulation that the holding 
condition only applies when s = mod N. As we vary 
N, we keep the amount of slack per stop, which we call 
<7, constant; the amount of slack per timepoint (that is, 
the amount of time a bus might actually have to wait 
while being held) is then No~. A computer simulation 
calculates the (dimensionless) buffer f3 of a "first bus" by 
determining the initial delay di for which the system 
moves from recovering (d^iooo < 10) and unrecoverable 
(di,iooo > 10: 10 and 1000 are both arbitrarily chosen). 

If we let B be the number of minutes corresponding to 
the dimensionless buffer j3, then according to Eq. [9] 
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FIG. 6: The amount of slack necessary per stop for a given 
amount of buffer B, both measured in minutes, as the number 
of stops between timepoints increases. The inset shows the 
amount of slack necessary per timepoint (that is, Na/B). 
Data points were found using computer simulation. 



Figure [6] shows how this ratio a / B depends on the pas- 
senger constant p! and the timepoint spacing N; for ex- 
ample, we find that for a busy bus route with N — 16 and 
[i! = 0.1, the ratio a/B — 0.21, so the route must have 
0.42 minutes' slack per stop, or a 6.7 minutes' slack per 
timepoint, to allow a bus to recover from a two-minute 
delay. For small passenger constants, the amount of slack 
per stop remains fairly constant, as if the holding from 
each stop can be "saved up" until the next timepoint. For 
larger passenger rates, however, holding becomes much 
less effective if the timepoints are spread too far apart. 



V. DISCUSSION 

When a single bus is delayed, it is able to recover if 
its delay is no larger than (3—1 (or B = a/ ji minutes), 
which we call the bus's buffer; the buffer is even larger 
for buses that trail already-delayed buses. 

If the number of passengers increases (such as during 
rush hour), the buffer will be reduced unless steps are 
taken to keep it steady, either by increasing slack, or by 



decreasing the unit boarding time (Eq. |2j; this is part of 
the reason that bus routes are generally more unreliable 
during peak hours. Since slack increases the time it takes 
to complete a route, one would like to find strategies that 
make it as small as possible: for example, it might be 
optimal to keep the slack proportional to the passenger 
arrival rate, keeping the buffer (3 = fi/a constant, so that 
busier stops get more slack. Alternatively, giving more 
slack to the stops immediately preceding or following the 
busier stops might be more advantageous; further study 
is required in this case. 

When it comes to choosing a holding strategy, both 
schedule and headway-based holding have advantages de- 
pending on the circumstances. When a series of buses ex- 
perience the same delay at a given stop, schedule-based 
holding can allow the buses to recover more quickly as 
long as the delay is no larger than the buffer. However, 
headway-based holding may be more appropriate in sce- 
narios with larger, random delays, as a very delayed bus 
will have a greater probability of trailing a slightly de- 
layed bus, and have a better opportunity to recover. 

One interesting, possibly counterintuitive result of our 
analysis is that buses on busier routes recover more 
quickly from small delays. The passenger constant seems 
to set the "timescale" of the route's behavior. 

Our next step will be to more fully investigate the effect 
of timepoints on these results. We have already shown 
that the necessary slack per stop increases dramatically 
when we only implement holding at every TVth stop. Fu- 
ture work will determine the exact nature of this relation- 
ship, and the effect that timepoints have on the behavior 
and buffers of subsequent buses. 

The model presented here is a "microscopic" model, 
dealing on the effects of delays at a single stop; how 
these single-stop delays interact with one another over 
the course of an entire bus route is also something worth 
investigating. Although engineers have spent a lot of 
time investigating this "macroscopic" regime, there may 
yet be basic truths which can be uncovered by a physi- 
cist's perspective. 

We thank Dr. Peter Furth of Northeastern University 
for useful conversations and his engineering perspective. 
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